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Abstract 

We have resummed all the (— boa s ) n contributions to the photon-meson tran- 
sition form factor F n ^. To do this, we have used the assumption of 'naive non- 
abelianization' (NNA). Within NNA, a series in (Nfa s ) n is interpreted as a series 
in (— boa s ) n by means of the restoration of the full first QCD /^-function coefficient 
— &o by hand. We have taken into account corrections to the leading order coefficient 
function and to the evolution of the distribution function. Due to conformal sym- 
metry constraints, it is possible to find the eigenfunctions of the evolution kernel. 
It turns out that the nondiagonal corrections are small, and neglecting them we 
obtained a representation for the distribution function with multiplicatively renor- 
malized moments. For a simple shape of the distribution function, which is close to 
the asymptotic shape, we find that the radiative correction decrease the LO by 30%, 
and the uncertainty in the resummation lies between 10% and 2% for Q 2 between 
2 and 10 GeV 2 . 
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1 Introduction 



We can expect that perturbative QCD (pQCD) works well to describe the process 
7*(gl)7(g2) — > ir° for accessible values of q\ < and q\ < 0. The form factor for this 
process is given by fjlj 



The leading term in a l/Q 2 expansion, being Q 2 the momentum transfer, is given by the 
integral in ([J). The dots, . . ., stand for higher twist contributions, which are sub leading 
in the l/Q 2 expansion. 

The coefficient function C(x, qi,q2), which accounts for the transition from photons to 
quarks, can be completely described within pQCD. It is known to one loop, and a detailed 
analysis can be found for example in [fj, |3], [|. 

The distribution function tp(x, fi 2 ) can be interpreted as the transition probability of a 
7T° with momentum P into two quarks with momenta xP and (1 — x)P respectively. Only 
its evolution with /x 2 is given by pQCD, but not its shape, or x dependence. Here, the 
situation is rather unclear, and there exist contradictory statements in the literature. Due 
to our present inability to extract it directly from experiment, we can only make some 
choice, (or guess) for ip(x,fi 2 ). The two most popular choices are the "asymptotic wave 
function", ip as = 6a; (1 — x), and the CZ model, tpcz = 30a; (1 —x)(l — 2a;) 2 . Very recently, 
new experimental data have appeared, and more are expected, f5|. It is expected that 
this will allow to constrain, and perhaps even extract with some accuracy, the distribution 
function (p(x,fi 2 ). In fact, one of the major goals of the study of the form factor F 7r7 is 
precisely to obtain more information on the distribution function. 

In order to extract as much information as possible from the experimental data, precise 
theoretical predictions for the coefficient function, beyond the one loop level, are needed, 
and it might even be necessary to include nonleading contributions in the l/Q 2 counting. 
In this work, we analyze the accuracy within which the leading order in l/Q 2 , that is, the 
integral in ([]]), can predict the form factor F nl within Naive Nonabelianization, NNA. 

The idea of NNA is based on the observation that corrections related to the evolution 
with the coupling can represent a source of potentially large perturbative coefficients. The 
extraction of these large contributions can give important information on higher order 
corrections. In QCD, this extraction can be done by evaluating the relevant feynman 
diagrams to leading order in the large Nf limit, and interpreting the series in (Nf<y s ) as a 
series in (— b a s ), restoring the full QCD /5-function coefficient — b by hand[]. Techniques 
to perform a summation of these large coefficients have been developed in 0. [7| . 

In this work, we will only calculate leading twist corrections, and use the ultraviolet 
dominance assumption || to estimate the higher twist corrections. This assumption is 
based on the observation that due to the (infrared) renormalon ambiguity, the leading 
twist result is affected by an intrinsic ambiguity, 

3 Here, for the (3- function, we adopt (3 = -b al + • • • , b — j^(^-N c - |iV/) 
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where A(y) is a calculable function that is completely fixed by the residues of the Borel 
integrand, see below, and y denotes the variables on which it depends. The ambiguity in 
@ is cancelled by another (ultraviolet) renormalon ambiguity in the higher twist contri- 
butions. According to the ultraviolet dominance assumption, not only the ambiguity, but 
the whole higher twist contribution is proportional to A(y). 

This work is organized as follows: In the next section, we briefly review the current 
state of art. In section 3, we compute the coefficient function within NNA. We find that 
the coefficient function has two 7i?-renormalon poles and that when one of the photons is 
on shell, there are additional ambiguities coming from the region x — > and x — *■ 1. These 
new ambiguities, which are related to the infrared region, lead to new power corrections 
to the form factor. In section 4, the NNA evolution kernel is presented. We obtain an 
expansion for the distribution function in series of Gegenbauer polynomials with the upper 
index shifted by bo<y s . We find that the nondiagonal part of the anomalous dimension 
matrix in this basis is much smaller than the diagonal part. In a first approximation, we 
can neglect these nondiagonal terms, obtaining multiplicatively renormalized moments. 
In section 4, we also present two models for the wave functions, and make some comments 
on the implications that conformal symmetry has. In section 5 we obtain the final result 
for the form factor in the NNA approximation. We present some numerical results for a 
simplified distribution function, where only the first term of the expansion is kept. In this 
case, the shape of the distribution function is close to the asymptotic one. In section 6 
we present our conclusions, and finally, we present two appendices with technical details 
of the calculations. 



2 The meson— photon transition form factor 

The transition from two photons to a tt° meson, 

7*(gi)+7fe)-vr°(P) (3) 

is described by the amplitude T 

T = e a ( qi )e^q 2 )4e apXp P x ^( qi - q 2 ) p F^{Q\ u). (4) 

Here, F nl (Q 2 ,uj) is the photon-meson form factor, e a (q\) and e' 3 ^) are the polarizations 
of the colliding photons, P = q\ + q 2 , —q\ > 0, — q\ > 0, Q 2 = —(qi — (fe) 2 /^ and u = 
(<7i — 4%) / (Qi + Q2) — 1 i s t ne parameter of asymmetry of the photons. If one of the photons 
(72) is real, q\ = and uj = 1. In experimentally accessible regions u is very close to 
one. 

In fig.|l] we have represented the process diagrammatically. The dominant contribution 
is given by figjlj.1, where a large virtual momentum flows through the subgraph containing 
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the two photon vertices. The other regimes correspond to a long distance propagation in 
the q2 channel ||. The second regime corresponds to the case where the large momentum 
flows through the central block containing a large virtual photon fig.|l[2. The third regime 
represents the situation when one of the quarks absorbs a large virtual momentum and 
carries almost all the momentum of the hadron and the second quark is soft, fig.|l|.3. Power 
counting predicts that the leading order for these contributions is 1/Q 4 . While for uo ~ 1 
all three regimes are important, only the first one is relevant for the situation in which 
both photons are off shell, to < 1. We will now discuss this contribution. 

For large Q 2 , the form factor F^Q 2 , to) can be expressed as the convolution of the 
coefficient function C(x, to, Q, fi) and the distribution function tp(x,fi 2 ) 



N 

~Q 2 



i j\r 
C(x,io,Q 2 ,/j 2 )tp(x,/j 2 )dx ee — C(x,to,Q 2 , fx 2 ) <g> <^(x,/i 2 ), (5) 

Q 



Here, N = /^(e 2 — e 2 ,) is a normalization factor, e g are the chargers of the quarks, 
/i 2 is the renormalization mass, or the scale that separates large from short distances, 
and f n is the pion decay constant, normalized to f n = 130 MeV, see for example fTO)] . 
The coefficient function can be calculated in perturbation theory from the hard parton 
subprocess 7* + 7* — > qq. At present, the coefficient function is known to leading twist to 
one loop accuracy, and in the limit to — 1 it has the simple form [B[ 



C(x,Q 2 ,/i 2 ) 



with x 



2x V 



In 



Ai 2 



(3 + 21nx) + lnx 2 - Ins- 9 
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x. The distribution function can be determined by the moments 



x}. 
(6) 



f n / x k tp(x,n 2 )dx 



(Pn) 



k+l 



(0\d l5 (n)(Dn) k u\P) 



(7) 



where is a light-like vector, n = 'y^n' 1 , the covariant derivative and \P) the 7r° meson 

1 

state with momentum P. By definition, tp(x,fi 2 ) is normalized to one, / ip(x,/j, 2 )dx = 1, 



and G-parity implies the relation tp{x) = ip(l — x). Its ji 2 dependence is described by the 
evolution equation 



d 



dfj, 2 

with some initial condition 



+ (3{a s 



_d_ 

da.. 



<p{x,n 2 



V(x,y\a s )Lp(y,n 2 )dy 



The kernel V(x, y\a s ) is calculable in pQCD, and can be expanded in series of a s : 



V(x,y\as) = ^V^(x,y)+^) 2 V^ 2 \x,y) 



+ 



(8) 



(9) 



(10) 
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Figure 1: The three regimes of that contribute to the process we are considering 
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It is known to two loop accuracy [IT], [JHJ. In the one loop approximation [IB 

1 A . 



VW(x,y) 



2C f 



x 



6{y>x)-\l- 

y \ y-x j 



[X 



x, y <r+y) 



[ID 



F(x,y)-S(x-y) [ F(t,y)dt. 
Jo 



Conformal symmetry, which at leading order manifests itself through 

y(l - y)vV)(x,y) = x(l - x)V^(y,x), 



(12) 



implies that the eigenfunctions that diagonalize the kernel (]TTJ) are Gegenbauer polyno- 
mials multiplied by xx. One therefore expands the distribution function in this basis: 



xx b n (n 2 

n=0 



<p(x,fi 2 ) 

The coefficients b n (fj, 2 ) are given by 

&n(/i 2 ) = b n {nl. 



2(3 + An) 3/2 



l + n)(l + 2n) 



2n 



2x) . 



(13) 



(14) 



Here 7^ are the eigenvalues of V^'(x,y), or the one loop anomalous dimensions of the 
multiplicatively renormalized operators: 



d^{h){nd)^C^{nD^/{nd) + )u 



(15) 



where (nd) + = nd + n d , 
eigenvalues are given by: 

7n("s) = 



1 n 



nD_ = nD — nD, see also the operators in (|7]). The 



O'.s- 



47T 

c F 



7« + 



1 2n+l -1 

+ 4V - 

;i + 2n)(l + n) p 2 j 



(16) 



and 6 = (H/3iV c - 2/3N f )/Air is the first coefficient of the QCD (3 - function. To 
next-to-leading order, conformal symmetry is broken by renormalization of the coupling 
and by renormalization scheme effects [|T4j], [IS|. The solution has a more complicated 



nondiagonal form, see for example ||, ||16|| . Assuming that the expansion in Gegenbauer 



polynomials (|13|) converges well, only the first harmonics are needed to obtain the model 
for the distribution function. The coefficients 6 n (A*o) i n (0) should be extracted from the 
initial condition <po(x), which is the low energy shape of the distribution function, and 
can not be calculated within pQCD. At present, there exist two popular models for <po(x). 
These are ip a s{x) = 6x(l — x), corresponding to the asymptotic distribution function in 
the leading logarithmic approximation, that is, only the first term in the expansion ( IB ) , 
and (pcz(x) = 30a;(l — x)(l — 2a;) 2 which has been proposed in |17| . In this model, the two 



first harmonics are taken into account. The second coefficient — 1 GeV ) = 3 has 
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been estimated using sum rules. In this model, the second coefficient is large and must 
be taken into account, while all the higher coefficients are assumed to be small and are 



neglected. In ||18|| ) the convergence has claimed to be slow, and therefore approximating 
the distribution function to the first few terms by ([13]) might not be justified. 

The relative contribution of all corrections depends on the choice of the model. In par- 
ticular, it has been suggested in |J that the discrepancy of the predictions for F nl (Q 2 ,u) 
within these models will be large enough to allow for an experimental discrimination. 



3 Calculation of the coefficient function 

To sum all (— boa s ) n contributions to the coefficient function, we have to calculate the 
coefficient function to leading order in the 1/Nf expansion, and perform the replacement 
2/3Nf — > 2/3Nf — ll/3iV c in the final result, according to the prescription of NNA. In 
our case, we have to calculate the one loop diagrams, but inserting a chain of n fermion 
bubbles. These diagrams yield factorialy growing contributions: 

r n ~ K(-a s b ) n n\n b (17) 



The convergence radius of the series in (jT7D is zero, and in order to perform a "summation" , 
Borel integral techniques are used. The bad asymptotic behavior of (|17|) will now manifest 
itself through renormalon poles in the integrand of the Borel integral, and a prescription 
has to be fixed to integrate over these renormalon poles. In this work we are going to 
make use of the principal value prescription. The result will depend on how we have 
integrated over the poles, being the ambiguity caused by the choice of a prescription 
known to be suppressed by powers of A 2 /Q 2 . This ambiguity is usually referred to as 
renormalon ambiguity, and it can be cancelled by taking the contributions of higher twist 
operators into account: if the renormalization of the lowest and higher twist contributions 
is performed consistently, an ambiguity free result can be obtained. 

In this work, we will evaluate the diagrams that contribute to the coefficient function to 
order 1/Nf, see fig.[§ The gluon line with a blob denotes the sum of all simple insertions 
of fermion bubbles, fig.^. The Born contribution, see fig^, is well known since long 
ago. We will use dimensional regularization to regularize the ultraviolet divergences of 
these diagrams. To obtain the renormalized contribution from each diagram, first the 
subdivergences, due to the fermion bubbles, and finally the overall divergence of the whole 
diagram, have to be subtracted. This we will do by following the technique presented in 
||. Details on the calculations, and the separate contributions of the diagrams can be 
found in Appendix A. Here, we simply present the final NNA results, together with some 
comments. We will present our results for the coefficient function in the following way: 

C(x,cu,Q7/i 2 ) = C (x,uj) + C (l-x,uj)+C 1 (x,uj,Q 2 /iJ 2 ) + C 1 (l~x,iu,Q 2 /fi 2 ) (18) 
Here Cq is the leading order Born contribution 

C {x,u>) = 1 —— (19) 

1 — UO + LXUJ 
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Figure 2: The diagrams that contribute to the coefficient function in the leading order 
approximation 
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_|_ v/V^^^A* _|_ . . . _|_ v/V^^^A/» ... y/V^^^/V _|_ 



Figure 3: The full blob represents the sum of all fermion bubble insertions 
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Figure 4: The two diagrams that contribute to the coefficient function in the Born 
approximation. 

In the next-to- leading order contribution, C\, a dependence on a s (/x) is understood: 
C F 



Ci(x,u;,QV^) 



-Co(x, to) x 



'0 



p.v. J e 
o 

l+a s 6o 



U 



'^c aM y ,,*«± u \ A -7(01*,.) 



(l-u)(2-u) 



+ 

7(u|x,a;) = 3 + 2 -Fi 
— 2 -Pi 



dA 



1-A 



[G(A|x,a;) -G(l|x,a;)] ^ , 



1,1 
2 + u 



2ux 



1,1 

2 +M 



1 + UJ 

2ujx 



2-u 



1 - uo 



2ujCq{x 1 uj) 1 + u 
2-u 2 



2ljCo(x, lu) 1 + u 



(20) 



(21) 



G{X\x,u) 



T(2A + 2) 



3T(2- A)r 3 (l + A) 



1,A 
1 + A 



2oox 



1-uj 



2 + 2A- A 2 
1 + A 
A 

2ujC (x, uj) 



+ 2^1 



+ 1 



1,A 
1 + A 



2ujx 



1 + u 



A 



2uCq(x, u) 



(22) 



where p.v. means that we use the principal value prescription to integrate over the IR- 
renormalon poles at u — 1 and u = 2. C parametrizes the renormalization scheme. In the 
MS scheme, C = ln(47r) — 7# + 5/3 [] (7$ is the Euler constant) and in the MS'-scheme, 
C = 5/3. The fact that we have to fix a prescription to integrate over the poles at u — 1 
and u = 2 induces an ambiguity 5Ci(x, u) in the coefficient function. It is well known 

4 Notice that our C, and the C used in [?|] have different globai signs. 



8 



that the ambiguity induced by a singularity at u = uq is power suppressed by (A 2 /^ 2 )" . 
In our case, for /i 2 = Q 2 , the ambiguity of the coefficient function will read 

Sd{x, u) = ±A 2 (x, uj)—± A 4 (x, u) I — ] (23) 



where A 2 (x, u) and A±(x,u) are the residues of the poles of the Borel integral, (|20|) , at 
u = 1 and u = 2 respectively. The /.R-renormalon ambiguity 5Ci(x, u) has to be canceled 
exactly by another ambiguity, the £7^/-renormalon ambiguity of the matrix elements of 
higher-twist operators. This means that higher orders in perturbation theory and higher 
twist contributions are inseparable. This fact can be used to obtain information on higher 
twist effects. According to the assumption of ultraviolet dominance, the full higher twist 
contributions are proportional to the UV— renormalon contributions, and the entire higher 
twist contribution can be included in the following way: 

N r 1 { A 2 A 4 1 

F^(Q 2 ,uj) =-^J o IC( Xi u,Q 2 ) + N 2 A 2 (x,uj)— + N 4 A 4 (x,uj)— + ■■■*> <p(x,Q 2 )dx 

(24) 

The dots denote higher power contributions. Within the assumption of ultraviolet domi- 
nance, the complete dependence on the kinematic variables x, u> is fixed by the calcu- 
lable functions A 2 (x, u), A 4 (x, u). The constants N t and their sign have to be fixed from 
experiment and it seems reasonable to expect their values to be of order one. In what 
follows, we set these constants to plus minus one and use this range as an estimate of the 
higher-twist effects. It should be kept in mind that this is a (perhaps raw) estimate. In 
Deep Inelastic Scattering, Ni ~ 2 [19], in support of our estimation. 

We now turn to the region that is accessible to experiment, that is, u 
only leading corrections, we obtain 

1 

2x ' 

RCV 1 



1 . Keeping 



Cn(xA) 



C x {x,l,Q 2 /^) 
C^Q 2 /^ 2 ) 



47T&0 



2x 



p.v. I I 





U 



e c n 2 



(25) 



l-«)(2-«) \2xQ 2 



l(u,x) -7(0,2;) 



+ 



1 



1 - A 



[G(X,x)-G(l,x)} 



(26) 



j(u,x) 



G(X,x) 



3+ 2 Fi 

T(2A + 2) 



" 1,1 






2 + u 


1 — X 


{ 







1 



x{2 -u)- 
2 + 2X-X 2 



1+u 



3T(2- A)r 3 (l + A) 1 1 + A 



+ (Ax - l) 2 F x 



' 1,A 






1+A 


1 — X 


) 







We have checked for u) — 1 that in the one loop limit these formulae are in agreement 
with (P). Consider Ci(x,Q 2 /[i 2 ). Using the simple identities 

1 111 



u(l-u)(2-u) 



1 

2u + 



(1-u) 2(2 -u) 



(27) 
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exp{-u/a s (n 2 )b }t u = exp{-u/a> s (/j?/t)bo} 



(2? 



we can rewrite the Borel integral in in the following way 



p.v. 



-u/asbo 







du 



it 



e c u 2 



[l-u){2-u) \2xQ 2 



l(u,x) - 7(0, x) 



(29) 



-u 



p.v . J exp 
o 

+7(0, x) lo| 



a s (2xQ 2 e~ c )b 
a s (2xQ 2 e~ c ) 



7(m,x) -l{0,x) 2 l{u,x) _ -f{u,x) 
u 1 — u (2 — u) 



du 



a s (u 2 



It is now easy to see that the integral diverges for x < e c A 2 /2Q 2 due to the Landau pole 
in the running coupling. This effect arises because our effective expansion parameter is 
A 2 /(xQ 2 ). This means that for small values of x, the entire power expansion in 1/Q 2 
needs to be resummed. A similar situation has been discussed recently in [pOj] . It has 
been shown that this resummation leads to new power corrections. It seems reasonable 
to assume that in our case these new power corrections are related to the fact that for 
to — 1, standard factorization breaks down for the higher twist contributions, and new 
regimes have to be taken into account, see discussion in Section 1. 



We will overcome this small x problem following the approach of f20 
convolute the coefficient function, 



We will first 

with the wave function, and afterwards do the 
Borel integral. The result will of course depend on the distribution function we have 
convoluted with, but in general, we will encounter a new singularity structure. For the 
asymptotic wave function, (p as (x) = 6x(l — x), for example, we will find singularities for 
all positive integers u, simple poles for u > 2, and double poles for u = 1 and u = 2. A 
very interesting situation arises with the NNA asymptotic distribution function, see next 
section for details. Here, the new singularities will be located at u = boa s (Q 2 ) + m, with 
m = 1,2,.... Despite the position of these new poles, the power corrections they induce 
are integer powers of A 2 /Q 2 : 



'A 



2 \ m+b a 3 (Q 2 ) 



'A 



2 \ m+l/log(Q 2 /A 2 ) 



(30) 



and our leading order result is in agreement with the predictions of the power counting 
rules. 



4 The distribution function 



4.1 The distribution function in the NNA approximation 

The pion distribution function is a phenomenological model function, and information 
about its shape should be taken either from experiment, or from nonperturbative calcu- 
lations. In perturbative QCD it is only possible to predict its evolution with u 2 using 
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the evolution equation (|S|). At the one loop level, conformal symmetry allows to find a 
basis of multiplicatively renormalized operators, flT5|). It follows that the eigenfunctions 
that diagonalize the evolution kernel, (|TT|) , are Gegenbauer polynomials multiplied by xx. 
This suggests to expand f(x) in terms of these eigenfunctions, (|13|). Since for n > 0, 
7^ > and 70 = 0, see (0), only the lowest harmonic in ( |13|) survives in the limit 
yU 2 — > 00. This leads to one of the most popular models for the distribution function, the 
asymptotic distribution function, <f as (x) = 6x(l — x). The two loop corrections break 
conformal invariance. The operators ( |T3"D get mixed, and the evolution kernel can no 
longer be diagonalized with the basis of Gegenbauer polynomials. As a consequence, the 
evolution of the distribution function is now much more complicated than in the previous 
case, see for example ||, [HJ. Below, we present the analysis of the distribution function 
within the NNA approximation. 

Our starting point is the lowest order NNA evolution equation: 



_d_ 



t, 2 9 



V[x,y\a s (ii 2 ))(p(y,n 2 )dy 



(31) 



The NNA evolution kernel reads 
V[x,y\a 8 (n 2 )] = V a (x,y) = 



T(2a + 2) 



Og_ a 

4?r F 3 r(2-a)r 3 (l + a) 



Q(y >x)[- 

\y 



Here, we have introduced the notation x = 1 



x 



a + 



1 



y- x y 



x, y 



y , a 



{x <-> x, y <-> y} 



1 + a s b and as usual 



[F(x,y)] + = F(x, y) — 5(x — y) Jq dtF(t,y). We have obtained this kernel combining 
the techniques presented in [^TJ and 0. 

The relevant diagrams are shown in fig.[5] Q This kernel, as in the one loop case, (|l2f) 
becomes symmetric after multiplication by (yy) a : 

(yy) a V a (x,y) = W(y,x) = W(x,y) . 

This fact allows us to obtain the eigenfunctions and eigenvalue of the kernel V a (x,y) f] . 



V a (x,y)<p n (y,n 2 



ip n {x,ii 



f 

Jo 



ip n (x,fi 2 ) 
A n (a s ) 

<~p n (p k dx 



-7 n (a s )v3 n (x,/i 2 ) , 
(xx) a ip n (x, fi 2 )A n (a s ) , 

C*2n + (1 — 2x) , 

r(l + 2a) (2n)! (1 



2a + An) 



r(a)r(l + a) (2a) 2n 



a + n 



(32) 

(33) 
(34) 

(35) 
(36) 



5 As a byproduct of our calculation we also obtained the evolution kernel P(Z) for the forward case. 
We agree with the result presented in pl| 

6 After finishing our calculations, Mikhailov published a work pif , where this possibility is also 
discussed. 
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1 



2 




Figure 5: The diagrams that contribute to the evolution kernel in the leading order 
approximation. For diagrams (2) and (3), the mirrors conjugate diagrams have to be 
added. 
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The eigenvalues ■y n (cn s ) read: 



7n(" s ) 



c f 



a 2 r(2a + 2) 



a(a + V 



4tt 3(1 + a) r(l + a) 3 r(2-a) 
2(1 + a) 



(a + 2n)(a + 1 + 2n) 



a: 



;i + a + 2n)-^(l + «)]>■ 



(37) 



Here, ^(2) is given by ^(2) = ^lnT(z). Notice that, as in the one loop case, quark 
current conservation implies 70 (ct s ) = 0. From fl33|) we see that in our case it is natural 
to expand the distribution function in Cf„ 2+Qs6 ° Gegenbauer polynomials: 

00 

^(x,/i 2 ) = (xx) 1+a ^Y,U^Mn(a s )C 3 2 i 2+asb0 (l-2x) , (38) 

n=0 

compare with ([TBI). Substituting fl3"8] ) in the evolution equation (|3"TD and using orthogo- 
nality of and <y9fc, we obtain the following equation for the moments 6fc(/z 2 ): 



2 9 L 2 



dfi 2 



_d_ 



h(n 2 ) = -h{^ 2 )lk{us) -^2C M (a s )b e (n 2 



(39) 



£=0 



where we have introduced the mixing matrix Cki{ot s ), which arises due to the fact that 
now the eigenfunctions depend on a s (/i 2 ). Technical details can be found in Appendix B. 
In fl39|), 60 is the first coefficient of the j3 function, and should not be confused with the 
moment &o (/•*)> where an explicit /1 dependence is indicated. We obtain a triangular system 
of linear differential equations for the coefficients b p (fi 2 ), p = 0, 1, . . . , k. Introducing the 
vector 

/ &o(/i 2 ) \ 
bi(» 2 ) 



B 



bki/J 2 



V ! 



(40) 



/ 



our equations can be written in the following matrix form: 

B(/x 2 ) = + C)B(^) 



2 g _ b 2 9 
djj? s da s 



(41) 



Here, jd is a diagonal matrix, built of the eigenvalues (|37f), and C is the triangular 
mixing matrix of ([39]). The first N elements of ([It]) can be obtained exactly for arbitrary 
N. One first notices that ( |39"D implies that &o(a* 2 ) is constant because Coo = 7o = (due 
to conservation of the axial current). For 61 (/i 2 ), (B9I) implies 



/'" 



_d_ 



- b al 



_d_ 



h (/i 2 



■[ji(ot s ) + Cn(a s )]6i(/i 2 ) - C 10 (a s )b (fi 2 



(42) 
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which can be solved using well known techniques. This procedure can now be repeated 
for &2O 2 ), h(fi 2 ), . . .. 

The general solution to (|39|) is given by 

B(^ 2 ) = U(n 2 ,fi)BQj*) (43) 

with the evolution matrix 

{V da' * ) 

U^\nl) = P a exp J ^Ma') + C(a')]\ (44) 

U(lilti) = i. (45) 

It follows that 

k 

h(v 2 ) = £ {^(/^o)L>G«o), (46) 

k'=0 

and we obtain the following expression for tp(x,fi 2 ): 



00 

^(x,/i 2 ) = (xx) 1+ ^£{^( M ^ / xg)} fefe; 6 fe K / x 2 )A fe K)^f +as6o (l-2x) , (47) 

fe=0 

where we sum for k' < k. 

We have found numerically that the matrix elements of the nondiagonal part of the 
mixing matrix C are much smaller than the diagonal matrix elements of the sum 7 + C 
for all resonable values of /1 2 . We plot the matrix elements for fi 2 = 5 GeV 2 in fig.^|a 
and fig.|B|.b. The diagonal elements are clearly bigger than the non diagonal ones. This 
happens for the whole range of relevant fi 2 . 

We can therefore solve equation (f|l]) by iterations with respect to the nondiagonal part 
C . At leading order, neglecting the nondiagonal part of the mixing matrix, all coefficients 
b k (fi 2 ) renormalize multiplicatively, and the distribution function has the following form: 

00 

y,(*,/i 2 ) = (xx) 1+b ^Y,U^)M^s)Cf +aM (l-2x) , (48) 

fc=0 

bkifi) = b k {n )exp\ — dx\, (49) 



where 



C kk {a s ) = (a s b ) 2 (\ + ol + 2k) - ip (a + ^) } . (50) 



4.2 Models for the distribution function 

We have already mentioned that the distribution functions is not predicted by perturbative 
QCD. In this subsection we discuss some models for the distribution function. For any 
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Figure 6: 

0, 1, 2, • • • , k — 1 and /i 2 = 5 GeV 2 . In ||.b we compare the nondiagonal elements, |6|.a, 
with the diagonal ones, 7 m + C mm , m < 15. The effect of the nondiagonal terms on 
the evolution of &i(/U 2 ), ^(/W 2 ) is displayed for (6 ,&i,^2)(l GeV 2 ) = (1,0,0) (§|.c), and for 
(b , 6i, 6 2 )(1 GeV 2 ) = (1, 1.7, 1.6) (|.d). 

physical distribution function, the first coefficient bo{fi 2 ) of the expansion (ESp is 1, as it 
follows from normalization 



1. 



(51) 



To extract the coefficients bk{$) for > 0, some information about the low energy shape 
of the distribution function is needed. At present, some sum rule estimates are available. 
The second moment 

[ (2x - 1) 2 ip(x, nl)dx = 0.35 {nl = 1 GeV 2 ) (52) 
J o 

has been estimated in [17], |T^j, and the following estimation for tp (1/2) can be found in 

ip(l/2, /i 2 ) = 1.2 ±0.3 (53) 
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For the second moment, the asymptotic distribution function predicts 



ip as (x)(2x - If 



0.2, 



(54) 



and for x = 1/2, its prediction is <p as (l/2) = 3/2. 



The estimation ( p2f) would imply &i(/Xq) = 11 in the one loop expansion fll3D , and 
^i(A*o) = 1-7 in the NNA expansion fl4*8|). The coefficient in the NNA expansion is 45% 
bigger than the 1 loop coefficient. To be consistent with the central value of fl5^), we need 
at least the third coefficient of the expansion of the distribution functions, that is, 



<p(x) = [x(l - x)] 1+b ° a ° (A (a s ) + MiC* 2 3/2+fe ° Qs (l - 2x) 



b 2 A 2 Cl /2+boas 



2x)) (55) 



in the NNA approximation, and ( j55|) with a s = in the one loop approximation. In ( 55 ) . 
Ai{a s ) are given by fl3"5"|). In the one loop expansion, (|T3]),this implies 62(^0) = an d m 
the NNA expansion, (|48|) ^2(^0) = 1-6- We have represented graphically both functions 
in fig.0. We agree with |1| that the higher harmonics remove the oscillations and make 
the function smoother. 




Figure 7: Several distribution functions. Solid line: asymptotic distribution function. 
Dot-dashed line: one-loop distribution function obtained by combining sum rule results 
(|52j) and the prediction (|53|). Dashed line: NNA distribution function, (p5|) obtained as 
above. Here, $ = 1 GeV 2 . 

The evolution of the moments bk{^ 2 ) with fi 2 is given by (P§p. It can be shown the 
in the limit \i — > 00, all &i>o(A i2 ) tend to zero. This follows from the fact that at the one 
loop level the NNA approximation reduces to the exact one loop approximation. Then 
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Figure 8: Solid line: asymptotic distribution function. Dashed line: Asymptotic NNA 
distribution function fl5T)|) for Q 2 = lOGeV 2 . Dot-dashed line: the same for Q 2 = lGeV 2 . 



only the lowest momentum, bo(fj, 2 ) will be relevant in the limit \i — > oo. This leads to the 
asymptotic wave function [} 

^ as{x) = J^ P(2^ °aj [x(1 ~ X)]1+b ° aa = 6X(1 ~ X) - (56) 

It is interesting to see the difference between the asymptotic shape and the lowest order 
NNA harmonic, which depends on a s . In fig.[| we can see that the lowest NNA harmonic 
is a little bit narrower than (p as (x). 

The effect of the nondiagonal terms Ck i on the evolution of the coefficients of the 
two models presented in this subsection has been displayed in fig.(||).c and fig.(||).d. In 
fig.(P).c we have solved the differential equation (|39|) exactly for the initial condition 
corresponding to (|56|), that is, &o(l GeV 2 ) = 1, &i>o(l GeV 2 ) = 0. Here, the diagonal 
solution ( [49D predicts &i(/i 2 ) = b 2 (fi 2 ) = 0, while the exact solution is negative, and quite 



small. In fig.(|5]).d we have used the initial conditions the nondiagonal elements i for 
the initial condition b (l GeV 2 ) = 1, &i(l GeV 2 ) = 1.7 and 6 2 (1 GeV 2 ) = 1.6, which 
is obtained by combining a sum rule prediction for the first moment, JTBL [17 ], and a 



7 Another way of showing that 6i>o(£t 2 ) tends to zero in the limit [i 2 — > oo is by induction on i: First, 
we solve (p2|), and show that bi(p?) x vanishes in the limit fi 2 — > oo. Then, the equation for 62 can be 
solved, implying that 62 also tends to zero. This is then repeated for 63, 64, • • • 
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prediction for (p(l/2), |Tj|. The moments 6fc(/i) 2 are the matrix elements of the operators 

J 7 5(^)(^)^C 2 3 / 2+60Os (n J D_/(^) + )n. (57) 

which renormalize multiplicatively. It might be convenient to redefine the one loop oper- 
ators ( |15|) in order to extract from the diagonal part all the main corrections. 



4.3 Conformal symmetry constraints 

We will now make some comments on the diagonalization of the evolution kernel and its 
eigenf unctions based on conformal invariance arguments. We start with the observation 
that at the one loop level we have conformal invariance, which gets broken when we go to 
higher orders. Conformal symmetry is broken by the introduction of the renormalization 
scale fi, and due to renormalization scheme effects [[15]. The breaking due to the renor- 



malization of the coupling are, to our accuracy, proportional to the first coefficient of the 
(3— function, while those due renormalization scheme effects are 1/Nj suppressed, and 
therefore lie beyond our accuracy. Despite the fact that conformal symmetry is broken, it 
still allows us to diagonalize the evolution kernel. This is equivalent to the diagonalization 
of the anomalous dimension matrix. The fact that the anomalous dimension matrix can 
be diagonalized at any order in perturbation theory follows from conformal invariance at 
the one loop level, where (|16| ) is implied. In general, the eigenvalues will be given by a se- 
ries in a s . Since there are no two eigenvalues with the same lowest order coefficient, fll6D , 
all eigenvalues are different, and we can therefore diagonalize the anomalous dimension 
matrix. 



Consider now the renormalization group equations for the operators defined in (|T5| ) 
for higher orders in perturbation theory: 

^-^°R = 1°R, (58) 

where Or are the renormalized operators, and 7 is the anomalous dimension matrix, which 
is only diagonal at the one loop level. Since it can be diagonalized, there is a matrix U(a s ) 
such that 

j(a s ) = U-\a a )^ D {a t )U{a a ) (59) 

where 70 is a diagonal matrix. Substituting this representation in fl58|) and redefining the 
operators 

R (a s ) = U{a s )0 R (60) 
it is easy to obtain the following equation for the new operators Or. 

M 2 ^A = {^(a.) - U(a s ) ^rV.lj I R (61) 

In the rhs, the coefficient in front of Or is its anomalous dimension. It is clear that to our 
accuracy, we can associate the diagonalization of the matrix 7 with the diagonalization 
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of the kernel V a (x,y). Then, the nondiagonal part, U(a s ) (/i 2 ^£7(as)J is the mixing 
matrix C(a s ) and the operator On(a s ) is given in (|5T|). 

In D = 4 — 2e dimensions, the /3-function has the following form: 

(3(e,a s ) = -ea s -b a 2 s -\ (62) 

and there is a critical value for the coupling (or fixed point ) a* such that 

/3{e,a* s ) = (63) 
To first order in the large Nf expansion, ( p2] ) implies 

3e 



2N f 



(Air) + 0(1 /Nf) (64) 



Consider now equation (^T|) at the critical point. Than nondiagonal part vanishes: 

U(a s ) L 2 -Lu-\ a .) J \ Qs=at = /3(e, a s )U(a s ) ^-Lfj- 1 ^ k=«: = 0, (65) 

due to (|63|), and 0#(a*) has a diagonal anomalous dimension matrix jr>(a*) at the critical 
point. For the operators (p7[) , taking into account that we have to perform the substitution 
2/3N f -> 2/3A^/ - ll/3iV c "by hand", and then a* = -e/b , we obtain 

R (a* s ) = d l5 (n)(nd)^Cli 2+boa HnDj(nd) + )u = d l5 (n)(nd)fC 3 2 i 2 - e (nD^/(nd) + )u 

(66) 

The last equation shows that the basis (|57| ) is conformal at the critical point because the 
classical conformal composite operator, built of two fermion fields ip, ip is given by [^]: 

^ 5 (h)(nd) 2 + n C^(nD„/(nd) + )ii , (67) 

where d^ is the canonical dimension of the fermion field. In D = 4 — 2e dimensions, 
= 3/2 — e. We have checked by explicit calculation that at the critical point, a* = 
2^(47r), and to leading order in the large Nf expansion, the operators (|66|) have a diagonal 

anomalous dimension matrix 7d(<^*)|. 

We conclude that conformal symmetry operators at the critical point and the operators 
On(a s ) that diagonalize the evolution kernel at leading order in 1/Nf are related: The 
operators On(a s ) can be obtained from the the conformal ones by the replacement (— e) — > 
b a s . We should notice that this result is rigorous at leading order in the 1/Nf expansion. 
However, we follow the prescription of NNA, and replace 2/3Nf — > 2/3Nf — 11/3N C . It 
seems reasonable to assume that the nondiagonal part of the anomalous dimension matrix 
induced by this will be small. 

8 Recently this anomalous dimension at the critical point has been calculated in |23] 
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5 The form factor 



Below we calculate the form factor (|) for uj — 1. In this case, (|) reads: 

F^(Q 2 ) =^f Q (Co(x) + Cx(x, Q 2 ) + C (l -x)+C 1 (l- x, Q 2 )) <p(x, Q 2 )dx , (68) 

where for simplicity we set n 2 = Q 2 . Co(x), Ci(x, Q 2 ) can be obtained from (^) and 
([2~6|), and (p(x,Q 2 ) is given by (|8|). 

We have already seen that the coefficient function (|26[) is ill defined in the limit x — > 
and x — > 1 because we expand in A 2 /(xQ 2 ) or A 2 /((l — x)Q 2 ). This induces uncertainties 
which are related to the nonperturbative structure of QCD in the infrared region and 
induce power corrections to the form factor. In addition to the power contributions 
from the 7i?-renormalon ambiguity, the form factor can get qualitatively different power 
corrections from the regions x — > and x — > 1. These new power corrections can be 
related to new contributions which correspond to the new regimes, see fig.|l}2 and [I]. 3. 
As discussed above, the leading contributions of these regimes are of order l/Q 4 and 
become essential in the limit u — > 1. To perform a phenomenological analysis of the 
power uncertainty related to the regions x — > and x — > 1, we first observe that the all 
the ambiguities come from the Borel integral of Ci(x, Q 2 ), see (p6|), and interchange the 
order of integration, that is, integrate first over x, and then over the Borel parameter u. 
Using the property ip(x) = ip(l — x) we rewrite (|68|) as: 



(C (x) + C^x.Q 2 )) <p(x,Q 2 )dx . 



(69) 



The complete integration over x can be done analytically. In principle, it is only necessary 
to interchange the order of integration in the Borel integral in (p6|), because the other 
integral is well defined. We obtained an expression for the form factor in which all 
ambiguities are related to the poles of the function in the Borel integral. The integration 
over x yields (below a = 1 + bo<y s ): 



Q 2 ^ {2a) 2 k (a + k) I 4tt6 



(70) 



M k (a s 



+ 



, r -u/boa a du I u (C-ln2) 2I> - U ) 

' J u \ (l-u)(2-u)r(a) 



9 k (a s ,u) - 9 k (a s ,0) 



1 - A 



{Pk{® s , A) - Pfc(a s , 1)} 



(71) 



9, (a 



x- 



,u 



,{l + u) 2k T(l + 2a + 2k) (l + u) 2k (2-u)(l + u)(a-u) 



(2k)\ 

T(l + 2a + 2k) 
T(2 + 2a - u + 2k) 
2T(l + 2a + 2k) 



T 1 + 2a 



u + 2k) § 



r(l + 2a - u + 2k) 

1, 1, 1 + a + 2k 
2 + u + 2k,2 + 2a 

(l + u)t 



(2k)\ 



(l+u + 2k) 



3 r 2 



U 



2k 



II 



2k 



[l +U + 



3-^2 



1, 1, 1 + a + 2k 
2 + u + £,l + 2a 



u + 2k 
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Pk{a s , A) 



+ 



T(2 + 2A) 

3r(2 - A)r 3 (i + A) \~ 

a\ 2 

(\ + 2k)(l + 2a + 2k) 



2 + 2A - A 2 



+ 



1 + A 

1, A, 1 + a + 2k 



3-^2 



1 + A + 2k, 2 + 2a + 2k 



2k 

E 



A 



A + 



3-^2 



1,A, 1 + a + 2k 

1 + A + L 1 + 2a + 2k 



(72) 



Where 3 F 2 [---|1] is the standard hypergeometric function (3,2), {z) n = r( ^^ is the 
Pochhammer symbol, and CV = {N 2 — 1)/(2N C ). We have absorbed all the singularities 
in the factor T(l + boa s — u)/[T(l + boa s )(l — u){2 — u)} of the Borel integral in (|7lD . 
There are two renormalon poles at u = 1 and u = 2, and an infinite number of new 
"small-x" poles at u = 1 + &o a s + wi, m = 0,l,---. We will integrate the poles following 
the principal value prescription. The final result will be prescription dependent, but this 
dependence is known to be power suppressed. 

To estimate the size of the higher power corrections, we will follow the recipe of Section 
2: A pole at u = uq induces a power correction of order (A 2 /Q 2 ) u °. The entire power 
correction to the form factor can written in the following way: 



AF^(Q 2 ) = ±7^(^ A ^)(^) +K ^ 2 ^{^-2) 
± ^ ({N 2 A 2 (a s ) + KM^e- 1 } (^) + 



2 \ l+b a s 



+ 



(73) 



In the last identity we have used the one-loop expression for the running coupling. By the 
dots, we denote the contributions of the higher poles, u = 2, and u — 1 + b a s +m, m = 
1,2, ■ • •. The contribution A^iA2(o;s) is related to the first renormalon pole, u — 1. The 
function A2(a s ) is build of the residues of the integrand of the Borel integral, and can 
be calculated, and N\ is an unknown number of order one. The second contribution is 
related to the "small-x" poles. Its structure is the same, K\ is an unknown constant and 
5i(a s ) a calculable function. We finally obtain: 



N /~ (2A;)! (l + 2a + 4k) 1 ,^ 2 , r „ x, 
Q \t^o ( 2a )2k {a + k) 

± {N 2 A 2 (a s ) + K^a^e- 1 } (^) +■ 



(74) 



Here, we are neglecting l/Q 4 and higher power contributions. To make a numerical 
estimation, we only take, for simplicity, the first term of the expansion of the distribution 
function fl56l) 



T{4 + 2b a s ] 
TV2 + b a s ) 



■[x(l — x)] 



l+6o«s 



a. 



= « S (Q 2 ). 



(75) 
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As discussed above, this approximation is very close to the asymptotic form. Then, in 
the expression for the form factor (|70D , we have to keep only the first term. This gives: 



F^JQ 2 ) = — - — — (l — — 7-M (<Xg)l ± power corrections 

-u/b a s ^_ j u(c-in2) T(a-u) 26 (a s ,u) . ' 

u f I» (l-u)(2-u) 0las,Uj 



N 

M (a s ) = ]>.r. j i 
dX 



(76) 



+ 



{p (a S} \) - po(a s , 1)} , 
1 1 — A 



9 (a s ,u) = r(l + 2a) 



(2-u)(a-u) 



r(l + 2a-u)(l + u) 



r(l + 2a-w) T(2 + 2a-M 

1,1, 1 + a . 
2 + it, 1 + 2a - u ' 



1,1, 1 + a 

2 + u,2 + 2a-u 



3-^2 



(77) 
(78) 

|1 

(79) 



Po(«s, A) 



At leading order, we obtain 



T(2 + 2A) 
3r(2- A)r 3 (l + A) 

1, A, 1 + a 
l + A,2 + 2a 



2 + 2A - A 2 



+ 



aA 



1 + A ' (l + 2a) 
1, A, 1 + a , 



OQ 2 



—3-^2 



1 + A, 1 + 2a 



T 7 y g 2 

To obtain the next-to-leading, or 1-loop result, we expand ([75|) in series of a s : 
Q 2 ) = Qx(l - x){\ + 6 a s (Q 2 )[5/3 + ln(x(l - x))}} + 0(a 2 s ). 



JO) 



(81) 



The coefficient function (|26| ) in the one loop limit gives the well known one-loop formula 
(PD, and to next-to-leading order in a s , we obtain: 



7 NLO t 

7T7 



AT r /•! 
« 2 ' = gil' + W. 1 



1 



+ ln[x(l — x) 



dx 



+ C, 



47T 7o 



;i-x) 
1 



ln2(3 + 21nx) + ln 2 (a;) - ln(l 



dx 



3^ (l " ~ 3 bo* s ~ BC F % 



52) 



The contribution of higher order radiative corrections in the NNA approximation is given 
by ([761). Numerical results are presented in fig.|9|. For a s (Q 2 ), we use a s {m 2 ) = 0.32, N j 



f 



4, and take into account that in the MS'-scheme C = 5/3. We see that the one loop 
correction decreases the leading order result by 20%. Higher order contributions, together 
with the one-loop correction, decrease the leading order by 30%. This estimation is in 
agreement with the estimations for radiative corrections suggested in 11. 
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Figure 9: The form factor F nl . The constant line (Q 2 /N)F nl = 3 corresponds to the 
lowest order prediction, ( jgQ|) . The next to leading result, (|82|) is plotted with long dashes, 
while the dots represent the NNA prediction, (|76|), and the two solid lines represent the 
uncertainties in the resummation. The dot-dashed line is the semisum of the NNA result 
and the next to leading prediction, and we expect the real result not to lie far away from 
it. 



The difference between our distribution function and the asymptotic one is very small 
and, for example, a discrimination in the description of the NLO correction to the form 
factor will be very difficult. Let us remind that in || it has been assumed that the 
corrections due to the evolution of the distribution function are so small that they can be 
neglected. 

To obtain the power corrections, we set the unknown constants K 2 and N 2 in (|73| ) 
equal to 1. We find that in our approach the power corrections give an ambiguity of 
order 10 — 2% for 2GeV 2 < Q 2 < lOGeV 2 . Below we give some numerical results in the 
following form: 



N- l Q 2 F^{Q 2 ) = F^ LO {Q 2 ) - (numberl) ± (number2), 



where "numberl" is the contribution of the higher order radiative corrections in the NNA 
approximation, and "number2" is the power correction (ambiguity in the summation of 
the perturbation series). 

For the last number, see (|73|), we set N% — N4 — K 2 = K4 = 1 (in this calculation 
we take into account small terms of order (A/Q) 4 ). Such choice fixes the relative sign of 



renormalon and small x ambiguities which is also unknown. In Fig.|10] we have plotted the 
separate contribution of the first 4 poles to the ambiguity. It is interesting to notice that 
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Q 2 (GeV) 2 




2 


2.20-0.36±0.19 


3 


2.27-0.36±0.14 


4 


2.31-0.35±0.12 


5 


2.34-0.34±0.10 


6 


2.36-0.33±0.08 


7 


2.38-0.32±0.07 


8 


2.40-0.31±0.07 


9 


2.41-0.30±0.06 


10 


2.42-0.29±0.06 



Table 1: Some numerical results for the form factor 



the first renormalon pole ambiguity, due to the u — 1 singularity, and the first small x 
ambiguity, due to the u = l+a s b have different sign, which results in a partial cancelation 
of the ambiguity. This cancelation suggests that the contributions of the matrix elements 
of the higher twist operators will partially cancel with contributions coming from photon 
emission at large distances. 

Other choice, such as for example N 2 = —K 2 = 1, which would lead to a bigger 
uncertainty (40% at Q 2 = 2 GeV 2 and 10% at Q 2 = 10 GeV 2 ), could also have been 
adopted. At present we have no rigorous arguments to choose the relative sign. Any of 
these definitions can only provide us with an estimation for the uncertainty. 

As a rule, NNA resummations give an exceeding estimation for the higher order contri- 
butions, see for example [0] , and the real curve for the form factor might lie a little above 
the one suggested by the NNA approximation in fig.^| but below the one loop approxi- 
mation. We have plotted an estimation for the real curve in the fig.|9] with long dashed 
(short gaps) line. 

To understand better the shape of the distribution function, we have to consider at 



least the model of the distribution function (55) which was obtained from the consideration 



of the sum rule constraints. This is work in progress, [E3 



6 Conclusions 

We have obtained the photon 7r°-meson form factor at leading order in 'naive nonabelian- 
ization' (NNA). We have computed the coefficient function and the evolution kernel that 
governs the evolution of the distribution function, in the NNA approximation. 

To evaluate the coefficient function and the evolution kernel, all leading contributions 
in (— bo<y s ) n have to be computed and resummed. For the coefficient function, Borel 
integral techniques have been used to perform the summation. Due to the presence 
of (IR) renormalons in the integrand, a prescription to perform the integration has to 
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Figure 10: The separate contribution of each pole to the uncertainty. 

be chosen. We have used the principal value prescription. The assumption of ultraviolet 
dominance of higher twist matrix element has allowed us to obtain an estimation of higher 
twist contributions (ambiguity in the resummation). We have found that when one of the 
photons is on shell, the coefficient function has new singularities. These new singularities 
are caused by new IR regions, that appear when one of the quarks that build the pion 
carries almost no momentum. The power correction induced by these new singularities 
have been taken into account. 

We have found that conformal constraints allow to find a basis of operators that di- 
agonalizes the evolution kernel in the NNA approximation, but the anomalous dimension 
matrix is nondiagonal. The nondiagonal terms of the anomalous dimension matrix have 
been found to be much smaller than the diagonal elements and we neglected by them. We 
have obtained an expansion of the distribution function in series of Gegenbauer polynomi- 
als with multiplicatively renormalized moments. This form is very convenient to analyze 
the effects of radiative corrections. Radiative corrections make the shape of the first term 
in the expansion of the distribution function a little bit narrower as compared to the one 
loop term but, we find that this effect is small. We have calculated the convolution of 
the coefficient function and the distribution function, and obtained the form factor as a 
sum with unknown nonperturbative coefficients, which depend on the shape of the wave 
function in the low energy region. In our numerical analysis, we only keep the first term 
in the expansion of the distribution function. We obtained that the radiative correction 
decrease the leading order by 30%. This is consistent with the estimation for the radia- 
tive correction by ||. The value of the power suppressed corrections (ambiguity in the 
resummation) depends on the relative sign of renormalon and small x ambiguities. It will 
be smaller for the case when they have different sign and predict the uncertainty for the 
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value of the form factor from 10% at Q 2 = 2GeV 2 to 2% at Q 2 = lOGeV 2 . The differ- 
ent sign of the ambiguities implies that nonperturbative effects coming from the matrix 
elements of the higher twist operator and soft photon emission is partially cancelled. In 
case of the similar sign the uncertainty for the value of the form factor is bigger: from 
40% at Q 2 = 2GeV 2 to 10% at Q 2 = lOGeV 2 . Clearly, an investigation of the influence 
of the higher harmonics in the expansion of the distribution function is needed. This will 
allow us to understand better the structure of the distribution function in the low energy 
region. 
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Appendix A 

All calculations have been performed in dimensional regularization in D = 4 — 2e dimen- 
sions. For the gluon propagator, we use the Landau gauge. To solve the 75 - problem, we 
have followed the prescription of 0, where an anticommuting 75 was only used in the box 
diagram. The unrenormalized contribution of each diagram Di to the coefficient function 
can be represented in the following way: 



A 



~£ E ( ) (n+1) a^(/x 2 )(-^[(n + l)6,6] , (A.l) 



n=0 



pr , (, /i 2 V ( r(4-2e) V" r(i-e)r(i+tt)r(i-„) , 

where n is the number of fermion bubbles in the gluon propagator, Q 2 = — (gi — <?2) 2 /4 > 
0, Pq = t^Nj /Ait is the fermion part of the full (3$ = —b , and only the factor Ri(e, A, x, u) 
is different for different diagrams. We will refer to the contribution of the g^ v term in the 
gluon propagator as the "gauge invariant" contribution, and to the contribution from the 
k^k v j k 2 part as the "gauge dependent" contribution. The total contribution of the "gauge 
dependent" cancels, and we will in what follows only deal with the "gauge invariant" ones. 
In what follows, the notation 

a = l-e, F(z) = 2^1(1, «; 1 + a + u\z) (A.2) 

where iF\ is the hypergeometric function (2,1), will be understood. 

The contribution of the crossed diagram of R will be denoted by R, Co(x,u>) denotes 
the Born contribution, see (|19|)). The diagrams are numbered according to fig.@: 

Rx+Rx = C 1+ "(x, u) 2C * 2 ^ + U) + {x <-> 1 - x} , (A.3) 

1 + a — u 
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R 2 + R 3 = Cl +U {x, u)l 2(1 -a)(a + u) + | a( ° + " ] 



+ a — u 



f 2au(a — u) , \ „ / 2l^(1 — x) \ ] , , 

+ (TT^- 2(a + " ) j F l^^)} + {a;wl " a:} 

fl 3 + fl 2 = Cl + '{x,w) {2(1 - a) (a + u) + ^ + 

[ 1 + Q! — U 

(2au{a — u) . ,\ „ /— 2uix\ 

+ — -i J --2(a + u" ] ' ' 

\ 1 + a-tt v 



From here, the sub divergences and the overall divergence have to be subtracted. In 
25|j a technique to perform this in the MS scheme was developed. We will now extend 
their method to arbitrary schemes. First, one observes that a general diagram with n 
fermion blob insertions into the gluon propagator can be expressed as ( |A.1| ): 



C<»> = (-D" + 'g^ (n ( +$„" +1 g([" + !]«.«) (A.4) 

where E([n + l]e, e) is a continuous function for e = 0. To subtract the sub divergences of 
the fermion bubbles, we also have to consider diagrams where at least one of the fermion 
bubbles is replaced by a counterterm: 

" C ^ £ U J (n + i-fcy*!-* ^ + 1 - fc]e ' e) T" J (A ' 5) 



In ( |A.5[ ), we have (*) diagrams (permutations) containing n — k fermion bubbles, (|A.4j) 



and counterterms. The sub divergences are cancelled by summing from to n. In the 
MS scheme, only the poles in e are subtracted, and Qms — 0, but in other schemes also 
finite parts are subtracted, and in what follows, we keep Q 7^ 0. The structure of ( |A.5| ) 
suggests to define 

E{[n + l]e, e) = e {n+1)n * E n ([n + l]e, e) (A.6) 
In terms of Eq, ( |A.5|) reads 



_ e ^ Cp ^ £ t A (-«sti) n - k mn + , _ k] e) ( «M\ {A7) 

47r few (n + 1 - k)e n+1 ~ k U i,J \ e J K ' 

Now, E([n + l]e, e) is expanded in the following way: 

00 00 
E n ([n + l]e,e) = J2E J (e)([n + lW, E (e) = s(e) = (A.8) 

j=0 i=0 



We now use ( |A.8|) to expand Eq([ti + 1 — k]e, e). Keeping only terms that do not vanish 



in the limit e — ■> 0, the contribution of the n fermion blob insertion, after subtraction of 
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the subdivergences, reads (see [^] for details): 

- e^)%^ + {-lYn\E n+l {e)e"A (A.9) 

Here, we have used 

r l/(n + l) J = 
]T(-l) n - fe (n + 1 - A;)- 7 " 1 = I 0<J<n + l (A.10) 

fc=o [ (-l) n n! J = n + 1. 

The second term in ( |A.9|) is already UV finite and no subtraction needs to be done, but 
the first one is UV divergent, and the poles in e, together with some finite pieces have to 
be subtracted. To subtract the finite pieces, a prescription has to be fixed: we will cancel 
those finite pieces that come from the expansion of the exponential. This is equivalent 
to taking a s e Qe , and not a s as an expansion parameter. The overall divergence can now 
easily be subtracted, and we find the UV finite result 

- C F ^-(aM) n (fjy + (-l) n n\E n+1 (0)) (ATI) 



In the first term in ( ATI ), n can easily be summed from to oo: 

r a 'f n Rf^f 9n+i\ C F f»tf g(\) - g(0) 
47r V n + 1 / 47T/35 Jo X 

To sum the second term, we use an integral representation for n\: 



dX (A.12) 



4:71 Jo ^ n AnBk Jo u 

(A.13) 



n=0 4tt^ Jo 

Applying this technique, we obtain: 



[D 



CU 



i R 



4tt/3 / 

a s f3 f 



du 



u 



Q 2 



2'j i (-u\x,uj) 
(l + u)(2 + u) 



ji(0\x,uj) 



(A.14) 



Q\gi(\\x,w) - gi(0\x,u)] 



where 



qAX\x,u) = lim EAu, X] = — — — -. rr 

y v 1 ' ; u-o * L ' J 6(1-A)r(l + A)r 3 (2-A) 



Ri(X,0\x,uj) , (A.15) 



ji(—u\x,u) 



{2 + u) 
2(1 -u 



■Ri(0, —u\x, u), 



28 



and C = 5/3 + ln(47r) — 7^ + Q parametrizes the renormalizations scheme. In the MS - 
scheme, C = 5/3 + ln(47r) — 7^ and in the MS - scheme, C = 5/3. Performing substitution 
/3q — > —b and rewriting the last integral in terms of — A|x, uj) = gi(\\x, u) we obtain 
our final result: 

+ / j-jlGiiXlx^-GiiOlx^)] 
1 



ji 2 e c 



Q 2 



2%(u\x,u;) 

;i-«)(2-«) 



(A.16) 




Appendix B 

Here we present some technical remarks related to the solution of the evolution equation: 



2 d 



b a 2 s 



d 

da. 



(f{x,fi 2 



V a (x,y)tp(y,n 2 )dy 



(B.l) 



Let us remind for convenience some formulae for the eigenvalues of the evolution kernel: 



V a (x,y)<f n (y,ii 2 ) 
<^ n (x,/i 2 ) 

V9 n (x,/i 2 ) 

A n (a s ) 

(fnfikdx 



-7 n (a s )v3„(x,/i 2 ) , 
(xx) a Lp n (x,fx 2 )A n (a s ) , 

G<2n + (1 ~~ 2x) , 

r(l + 2a) (2n)! (1 + 2a + 4n) 
r(a)r(l + a) (2a) 2n 

Ofcn j 



(B.2) 



a + n 



(B.3) 



where a = 1 + 6oa s . It is natural to expand the solution in a series of eigenfunctions ( B.2|) 
of the evolution kernel: 



tp(x,/j, 2 



(xx) 1+asbo 



Y.U^)M*s)cll 2+asb \i-2x) 



(B.4) 



n=0 



Substituting ( |B.4| ) in ( |B.1| ) and using orthogonality of the eigenfunctions we obtain the 
evolution equation for the moments &fc(/i 2 ): 



D RG b k (^ 2 ) + h n{^ 2 ) Wk ® D RG ip n } = -7 fc (a)6 fc (/i 2 ). 

n=0 



;b.5) 



where we have introduced the notation 



D 



d 



RG = H> 



b a 2 



d 



d\i 2 s da s 

Vk®^n= [ Vn^kdx 
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The second term in the lhs of ( |B.5| ) arises due to the dependence of <p n on a s (/i 2 ). We 
rewrite it in the following way: 

oo oo oo 

Kin 2 ) Wk ® D RG Lp n } = b n (fi 2 )D RG {if k g> ipn} - bn{^ 2 ) {D RG ip k } ® (fi n (B.6) 



n=0 



n=0 



n=0 



The first term in the rhs of ( |B.6|) vanishes due to ( B.3 ). Consider the second one. 

d_ c s/2 +boas {1 _ 2x) = _ k 

£=0 



Di«;y!< ■'■■"-{[ - 2x) = -Y,Cu{a s )Cl /2+b ° aa {l - 2s) . (B.7) 



In order to calculate the coefficients Cm, it is convenient to use the following formulae for 
the Gegenbauer polynomials PEfl : 

c» k ( t ) = { - )2{k ~ P)HPT{2k ~ p + u) ■ 



(2t) 



2p 



p=0 
P 



(2[k-p))\ P \ r(iz) 



(2p)! 

We have obtained: 

k 

Ckk>{a) = (b a s ) 2 Y 

p=k 



^ rv T(u)(2p-2k' + u) 

h ^- 2k,[t) T{k' + l)T{2 P -k' + v + l) 



k'=0 



(-) k ^T(l + b a s + k + p) (2k' + b a s + l) 

, (k - p)\(p - k')\ r(| + b a s + k' + p) (p + k' + b a s + §) 



x |^ f- + &o" s + + p^ - (b a s + 



In particular: 



C, 



(6 a s ) 2 jv> (2 + & o" s + 2A;) - tp (b a s + ^ 



Where ip(z) = £\nT(z). Substituting in and using (|B.3|) we obtain 



2 d L 2 



6 fc (p 2 ) = -& fc (p 2 )7/k(a s ) -5^CM(o! a )6/(^ 2 

£=0 



(B.8) 
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